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Gigamasers: the key to the dust-obscured star- format ion 
history of the Universe? 
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ABSTRACT 

We discuss the possibility of using OH and H2O gigamasers to trace the redshift distri- 
bution of luminous, dust-obscured, star- forming galaxies. It has long been thought that 
ultraluminous, interacting galaxies should host gigamasers due to their vast pumping 
infrared (IR) luminosity, the large column density of molecules available to populate 
the maser states and the turbulent motion of the gas in these dynamically complex 
systems which allows unsaturated maser emission. OH masers may thus be well-suited 
to the redshift-blind detection of ultraluminous and hyperluminous infrared galaxies 
(-^FiR > 10^^ L0) such as those uncovered by the SCUBA submillimetre camera. The 
bandwidth requirement is low, < 1 GHz for z = 1-10 (lower still if additional redshift 
constraints are available) and the dual-line 1665-/1667-MHz OH spectral signature 
can act as a check on the reality of detections. 
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1 INTRODUCTION 



The discovery ot a distant population ot luminous sub- 



successful, resulting in four r edshifts from a sam ple of fi fteen 
weakly lensed galaxies in the Small et al. (199^) samp le (Ivi- 



miUimeU'e (ijubmm) galaxios has rtivoluLionisti d our tuiTter 



mpii 

son et al. 1998|; Soucail et al. 1999|; Ivison et aT. 200C| , Prayer 



standing of galaxy formation and evolution ( Small et al 



199c 



1997; Hughes et al. 1998|; Barger et al. 1998|; Bales et al 



CarlUi et al. 2000D . Debate continues regarding the 



et al.. In prep), three of these have so far been confirmed 
as massive and gas rich systems through CO line mapping 
([Prayer et al. 1998|, 1999|, Knelb et al., in prep). 



relative Importance of obscured and unobscured star for- 
mation, the fraction of active galactic nuclei (AGN) In the 
submm galaxy population and the relationship. If an y, be- 
tween Lyman-break galaxies and submm galaxies (Adel 



However, the realisation that the majority of the submm 
popu lation have no pla usible optical counterparts (/ > 26, 
e.g., Small et al. 200C), has meant that the possibility of 



tween layman- breal; galaxies and submm galaxies (|Adel- 
berger fc Stcldcl 200(|: Chapman et al. 200C|; Peacock ct al 



200q; Bales et al. 2000|; Small et al. 2002| ); nevertheless. It 
clear ,that rest-frame far-iufrarpH (far-TR) energy that h 



a complete optical spectroscop ic survey has had to be dis- 



missed (e.g., Barger et al. 1999). Bven in the IR, only around 



half of the galaxies are identified by K < 22, often as ex- 
tremely red objects (BROs , I — K > 5 



Small et al. 199^ 



been re^^rocesseH hv Hnst and reHshifteH into the^wavehanH Gear et al. 2000|; Ivlson et al. 2000D, giving little hope to IR 



accessible to SCUBA (Holland et al. 1999) traces a galaxy 
population which makes a significant, and possibly domi- 
nant, contribution to the star-formation density at z > 1. 

It was realised from the outset that the most crucial 
piece of information required to derive the history of ob- 
scured star formation from the submm pop ulation Is its 
redshift distribution, N{z) (Blaln et al. 1999). Knowledge 
of N{z) breaks degeneracies In the models and allows the 
nature of the galaxies to be explored, most Importantly pro- 
viding estimates of their masses via observations of CO. Ini- 
tial efforts to target the first few SCUBA galaxies were very 



spectroscoplsts either. 

As in other branches of extragalactlc research, atten- 
tion has therefore been diverted to broadband photometric 
redshift techniques. Unfortunately, the classical analysis of 
optlcal/IR photometry has been shown to be misleading for 
even the few visible examples of these very dusty galaxies, 
due to the complex effects of dust on the optical spectral en- 
ergy distribution, a problem which Is further compounded by 
the low photo metric precision ava ilable for these faint galax- 
ies. However, CarlUl & Yun (1999) made a significant break- 
through with the realisation that one of the strongest cor- 
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relatio ns in observational astronomy - between far-IR and relationship being closer to Loh oc Lpnf, which Diamond 



radio ei lission ([Gondon 1992 ) ~ could be exploited to give an 
indication of redshift, based on the submm fluxes and sensi- 
tive radio detections of the gala xies. Several variants of this 
techniq ue have been developed ( Carilli fc Yun 200C|; Barg"er 



et al. (1999) have proposed to result from the admixture of 



unsaturated and saturated emission. 

Regardless of the precise charac ter of any I/qh/I/fir re- 
lationship, however, it is clear (e.g., Baan 1989|; Baan et al 



et al. 2300; Dunne et al. 2000), though, based on compar- 



1992b; Briggs 1998; Darling fc Giovanelh 200C) that pow- 



ison against the few submm galaxies with known redshifts, 
none are demonstrably better than the original. The tech- 
nique has been used to demonstrate convincingly that the 
median redshift of submm galaxies must be in the range 
2 < 2 < 3.5, with no significant low-redshift tail (Small et 
al. 2000; cf. Small et al. 1998; Lilly et al. 1999). 



erful OH masing is relatively common within the ultralu- 
minous IR galaxy (ULIRG) population, with upwards of 50 
per cent of those galaxies with Lfir > 10^^~^^ L© supporting 
megamasers (Loh^IOL©) or gigamasers (Lqh 10^ L©). 
Based on the assumption that the luminous sub mm galaxies 



Hoiv else can progress be made towards the determi- 
nation ot the N{z) tor the submm populationY Some have 
proposed that searches for CO rotational lines cen tred on 
subr nm galaxy po sitions are the way forward (e.g., Hughes 
200C; E lain 200C). The CO lines are certainly expected to 
be luminous, but the bandwidth requirements (~ 100 GHz) 
of such an approach are several orders of magnitude beyond 
the capabilities of current i nstrumentati on. Moreover, prior 
to the advent of ALMA (Brown 1999), observers will be 
reliant on current/planned 10-50m mm/submm single-dish 
facilities. As it has long been known that even inter feromet- 
ric detections of CO (where one benefits from stable base- 



are t he high- redshift counterparts of ULIRGs (Small et al 



2002), it is therefore quite reasonable to expect strong OH 



maser emission from a large proportion of the submm pop- 
ulation, which, if detectable, will allow the accurate deter- 
mination of their N{z). 

Although 22.235-GHz H2O megamasers show some- 
what different characterist ics to the OH maser systems (e.g., 
Braatz et al. 1994, 1996), similar arguments may be de- 



ployed in favour of their potential as redshift markers for 
luminous submm galaxies. H2O masers are more typically 
associated with AGN and are less useful tracers of star for- 
mation than their OH counterparts, but a significant nu mber 



of submm galaxies are known to harbour active nuclei (Ivl 



lines) i 1 galaxies with known redshifts are extraordinarily son ct al. 1998, 2000, Frayer et al., in prep). If detected, high 



diflicult ( JFrayer et al. 1998| j; there are obviously substantial resoluti on studies may yield information on their black hole 



technical difficulties to overcome when undertaking redshift- 
blind, single-dish CO searches. It is clear therefore that we 
should explore other avenues to measure the redshifts of the 
submm galaxy population and so understand more about 
their detailed properties. 

In this paper we propose an alternative method for de- 
termining the redshifts of submm galaxies, based on the ex- 



pectation that these dusty, ultraluminous galaxies will ex- 
hibit similar behaviour of their H2O and OH maser activity 
to that observed in luminous IR galaxies in the local Uni- 
verse. We first discuss the background to this proposal (§2), 
based on the properties of H2O and OH masers in the local 
Universe, before investigating, in §3, the feasibility of these 
observations using current and future instrumentation. Fi- 
nally, in §4, we state the main conclusions of this study. 

2 H2O AND OH MEGAMASERS 

Our proposed technique for determining the N{z) of the 
submm population is based on the expectation that lu- 
minous submm galaxies (typically ~ lO^^/i^g^L©, where 
/i65 = -ffo/65 kms^^ Mpc"^) exhibit similar scaling of their 
H2O and OH maser activity to that observed in luminous 
IR galaxies in the local Universe. Observations of OH emis- 
sion in local luminous IR galaxies, primarily due to the 
1665-/1667-MHz ground rotational state transitions, indi- 
cate a strong Loh oc ipiR correlati on between (isotropic) 
maser and far-IR luminosities (e.g., Baan 1989). This rela- 



masses ( Miyoshi et al. 1997 ), in addition to their redshifts. 
The bandwidth requirement for detection of H2O masers is 
only ~ 9GHz (vohs = 2-11 GHz for z — 1-10), similar to 
the bandwith planned for the new Very Large Array (VLA) 
correlator. In the case of OH masers this requirement drops 
below 1 GHz {uobs = 165-835 MHz for z = 1-10). Since there 
is no signifi cant submm galaxy population at z < 1 (3mai] 



et al. 2000), the techniq ue avoids con tamination by Hi emis- 



sion from local galaxies (Briggs 1996 ) and we can thus search 
a relatively clean part of the electromagnetic spectrum using 
existing interferometers with superb instrumentation, high 
aperture efficiencies and large collecting areas (e.g. VLA, 
Westerbork Synthesis Radio Telescope - WSRT, Australia 
Telescope Compact Array - ATCA). 

Inevitably, the difficulties with this new tech- 
nique will lie in its stringent sensitivity and dynamic- 
range requirements, as well as radio- frequency inter- 
fcrcnce in the 165-835-MHz (UHF) band (see, e.g., 
tittp : //www. atnf . csiro . au/SKA/intmit). In the following 



tionship has been ex plained in terms of a model first pro- 
posed by Baan (1985), whereby an OH population inversion. 



efficiently pumped by the far-IR flux, provides unsaturated 



amplificat ion of the background radio continuum; the Con- 
don (199"^ ) far-IR/radio luminosity correlation then leads to 



the observed q uadratic dependen ce of Lqh on I/fir. More 
recent work by Kandalian (1996) has claimed that the cor- 
relation is merely a consequence of Malmquist bias, the true 



section, we estimate the sensitivity requirements by extrap- 
olating from previous observations of H2O/OH megamasers 
and gigamasers in luminous IR galaxies. 



3 SENSITIVITY REQUIREMENTS 

Although the strength of extragalactic maser sources is of- 
ten reported in terms of the integrated line flux, we prefer 
to conduct our estimations using the peak line flux density 
as it is this quantity which determines whether a source 
is detectable with a given instrument. For a maser at red- 
shift z, emitting with a peak rest-frame isotropic luminosity 
density of L^/, the observed flux density Si, at a frequency 
v = v' /{I + z) will be given by 



= (1 + z) 



(1) 
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Figure 1. Contour map of logj^Q S,^, the peak observed flux den- 
sity, as a function of redshift, z, and peak rest-frame luminosity 
density, L^i . Also shown are observations of H2O (filled triangles) 
and OH (open squares) maser sources in luminous IR galaxies; 
the data for t hese Doints have been calculated from observation s 
pubUshed bv |Baan ct al. fl992a[b|). ^tavclv-Smith ct a (1992| ), 
Braatz et al. (1996| ) and |Darling fc Giovanelli (200o|, 200l| ). Where 
more than one observation has been made of a given source, the 
vertical lines indicate the spread of L^i about its average value. 



erwise unchanged, the observed flux densities would be 
~ 0.089 mJy, ~ 0.093 mJy and ~ 0.40 mJy, re- 
spectively. The paucity of points in Fig. ^ with Si, < ImJy 
illustrates that the detection of such faint sources is 
probably beyond the capabilities of current technology. 

However, the OH sources shown in Fig. |l] are all em- 
bedded in /7?^5-selected ULIRGs, and hence constitute a 
f ar-IR flux-limited sample, with a redshif t cutoff at 2: ~ 0.4 
( [Clemen ts, Saunders fc McMahon 199£ ). The population 
(e.g., |Rowan^Robmsoir"2000[ ) of hyperluminous IR galax- 
ies (HLIRGs), with I/pir > 10^'^ hg^ Lq, suggests that more 
powerful OH masers may lie undetected at redshifts z > 0.4. 
The HLIRGs are especially promising candidate hosts for 
OH masing: observ ations indicate that the y are powered by 
starburst activity ( Rowan- Robinson 2000 ) which may pro- 
yide the turbulence required fo r unsaturated masing to occur 
(Burdyuzha & Komberg 1990). Recalling that such unsatu- 



rated emission exhibits a quadratic Loh oc -LpiR behaviour, 
it is therefore possible that OH masers in HLIRGs may ex- 
ist with peak luminosity densities L„i approaching two or- 
ders of magnitude greater than the values shown in Fig. ^ 
for the ULIRGs. These immense luminosities would render 
putative HLIRG gigamasers detectable out to z > 4 at the 
5*^ ~ 1 mJy level, close to the sensitivities of present-day 
instrumentation, and within the grasp of facilities such as e- 
VLA, e- MERLIN and — ultimately — the Square Kilometre 
Array (SKA); however, issues relating to dynamic range and 
interference will need to be addressed. 



where Dl{z) is the luminosity distance ( Hogg 1999| ). The 
(1 -I- z) factor in this expression accounts for the linewidth 
narrowing due to the redshift, and partially offsets the 
quadratic drop-off in Sv with distance. 

In Fig. 0, we plot a contour map of logjQ Su as a func- 
tion of 2 and L^f , where we have adopted /les = 1, 0.m = 0.3 
and r^A = 0.7 for the evaluation of Dl{z). Overplotted in 
the diagram are the loci of a selection of the H2O and OH 
maser sources observed in luminous IR galaxies (see caption 
for references). Evidently, the H2O sources appear clustered 
at lower redshifts {z < 0.3) and lower luminosity densities 
(L^/ < 5 X 10^^ WMHz"^); the OH masers tend to be found 
at higher redshifts and luminosities. Whether this distribu- 
tion is intrinsic or due to selection effects remains unclear, as 
does the apparent correlation between z and L^i , although 
this is most likely a manifestation of Malmquist bias. 

The three OH sources in the diagram with 
largest L^,, correspond to /7?y4S'20100-4156 



9.0 
2.1 X 



X 10^"WMHz- 
10^" WMHz"^) 



2.0 X 10^°WMHz^ 



ffi^S' 120324-0707 
and IRAS U070+0525 
The latter, discovered by 
Baan et al. (1992a), is the most luminous gigamaser 
system currently known (inferred Loh ~ 1-05 x 10^ L©); 
however, this energy is distributed over a velocity width 
of ~ 2,400kms-\ which explains why ffi^S' 14070-1-0525 
exhibits a smaller L^, than the less lu niinous yet narrower- 



lined gigamasers in /j^Ag 12032 +0707 ([Darling fc GiovaneUi 
2OO1I2 



nd /7?^S'20100-4156 ( ^tavely-Smith et al. 19891 ). 
If these three sources were at redshift 2 ~ 3 (a typical 
value anticipated for the submm galaxies) and left oth- 



4 CONCLUSIONS 

OH and H2O megamasers are common constituents of the 
most luminous IR galaxies in the local Universe. The strong 
evolution in the population of dusty starburs t galaxies re- 
vealed by recent submm observations (e.g.. Small et al 



1997) should thus result in a population of distant galax- 
ies - submm-selected galaxies ('SCUBA galaxies') - hosting 
extremely luminous masers. These lines should be bright 
enough to be at the limit of detectability with current in- 
struments, but within the reach of e-VLA, e-MERLIN and 
ultimately SKA. 

We propose that the redshifts of submm-selected galax- 
ies, largely beyond the reach of optical and IR spectro- 
scopists, can be determined using interferometric searches 
for these maser lines. Maser searches have several clear ad- 
vantages over other methods: 

• the bandwidth requirement is small, < 1 GHz for 2 = 1- 
10 for OH masers, smaller still if additional redshift con- 
straints are availa ble (from their radio -submm spectral in- 
dices, for example [Carilli fc Yun 1999 ); 

• the instantaneous survey area is limited by the primary 
beam of the interferometer - several degrees for an OH line 
search with e-VLA, for example; 

• interferometry permits some rejection of local radio- 
frequency interference; 

• the position of an emission line can be pinpointed ac- 
curately within the primary beam, tying an emission line to 
a submm galaxy unequivocally; 

• the dual- line 1665-/ 1667-MHz OH spectral signature 
can act as an important check on the line identification and 
the reality of detections. 
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Armed with accurate redshifts for a significant propor- 
tion of the submm galaxy population we could test the pro- 
posal that SCUBA galaxies represent the massive progen- 
itors of present-day ellipticals, using measurement of their 
gas masses and fractions from interferometric CO observa- 
tions. The redshift distribution for the SCUBA galaxies de- 
rived from observations of megamasers would also remove 
the final ambiguities in interpreting the contribution of this 
population to the total star- formation density at redshifts 
of 2 - 1-5 ( iBlain et al. 1999| ). 
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